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The well-known analogy between a special limit of General Relativity and electromagnetism is 
explored in the context of the Lorentz- violating Standard-Model Extension (SME). An analogy is 
developed for the minimal SME that connects a limit of the CPT-even component of the electro- 
magnetic sector to the gravitational sector. We show that components of the post-newtonian metric 
can be directly obtained from solutions to the electromagnetic sector. The method is illustrated 
with specific examples including static and rotating sources. Some unconventional effects that arise 
for Lorentz-violating electrostatics and magnetostatics have an analog in Lorentz-violating post- 
newtonian gravity. In particular, we show that even for static sources, gravitomagnetic fields arise 
in the presence of Lorentz violation. 



I. INTRODUCTION 

In its full generality, General Relativity (GR) is a 
highly nonlinear theory that bears little resemblance to 
classical Maxwell electrodynamics. Nonetheless, it has 
long been known that when gravitational fields are weak, 
and matter is slow moving, analogs of the electric and 
magnetic fields arise for gravity [1[. These fields are 
sourced by a scalar density and vector current density, 
just as in electrostatics and magnetostatics. Further- 
more, in the geodesic equation for a test body, terms 
of no more than linear order in the velocity resemble the 
classical Lorentz force law arising from effective gravi- 
toelectric and gravitomagnetic fields [2|. Also, a well- 
known analogy exists between the precession of classical 
spin in a gravitational field and the precession of the spin 
of a charged particle in an electromagnetic field [3, |3| • 

In this work we investigate the fate of the standard 
connection between stationary solutions of the Einstein 
and Maxwell theories when violations of local Lorentz 
symmetry are introduced. Recent interest in Lorentz vi- 
olation has been motivated by the possibility of uncov- 
ering experimental signatures from an underlying unified 
theory at the Planck scale [a-0]- We examine the mod- 
ified Einstein and Maxwell equations provided by the 
action-based Standard-Model Extension (SME) frame- 
work, which allows for generic Lorentz violation for both 
gravity and electromagnetism, among other forces |8l-[l0| . 

In the so-called minimal SME case, the electro- 
magnetic sector contains 23 observable coefficients for 
Lorentz violation organized into two parts: 4 CPT- 
odd coefiicients in [kAp)^ with dimensions of mass and 
19 CPT-even coefficients in the dimensionless {kp)^^'^^ 
[3, [ill . The former set has been stringently constrained 
by astrophysical observations at the level of 10""*^ GeV 
p^ [l3| . The latter set has been explored over the last 
decade using astrophysical observations [iJI and sensitive 
laboratory experiments including resonant-cavity tests 



|15| . among others [la].^ Currently, constraints on these 
19 coefiicients are at the level of IQ-^'^ to 10^^^ ^ 

In the minimal SME gravitational sector, there are 20 
coefficients for Lorentz violation organized into a scalar 
M, two-tensor s'"^, and four-tensor t^^"'^^ [l0|. Within the 
assumption of spontaneous Lorentz-symmetry breaking, 
the dominant effects for weak-field gravity are controlled 
by the subset called s'"^ [lS|. These gravity coefficients 
have been explored so far in lunar laser ranging |19| and 
atom interferometry |20l - [23 |. while possibilities exist for 
other tests [2j], including time-delay and Doppler tests 

Since the CPT-even portion of the electromagnetic sec- 
tor of the minimal SME has 19 coefficients and the grav- 
itational sector, apart from an unobservable scaling m, 
also has 19 coefficients, one might expect a correspon- 
dence between the two sectors - an extension of the con- 
ventional analogy. Indeed, as we show in this work, there 
is a correspondence under certain restrictions. Thus it 
turns out that, under certain circumstances, Lorentz vi- 
olation affects classical electromagnetic systems in flat 
spacetime in a similar manner as gravitational systems 
are affected by Lorentz violation in the weak-field limit 
of gravity. As a consequence, some of the unusual effects 
that occur for Lorentz-violating electromagnetism have 
an analog in the gravitational case. In addition, from 
a practical perspective, it is quite useful to be able to 
translate analytical results in one sector directly into the 
other, as we illustrate toward the end of this work. 

We begin in Sec. |TT] by reviewing the basic field equa- 
tions for the gravitational and electromagnetic sectors of 
the SME. Next we explore the solutions to these equa- 
tions and establish the analogy between the two sectors 
in both the conventional case and the Lorentz-violating 
case in Sec. IIIII In Sec. IIV[ we explore test body mo- 
tion for both sectors and establish the connection in the 
conventional and Lorentz-violating cases. We conclude 
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^ For a thorough list of experiments, see the collected data tables 
inRef. [ij]. 



this work in Sec. |V] by illustrating the results with the 
examples of a pointlike source and a rotating spherical 
source, and we discuss some experimental applications of 
the results. Finally in Sec. IVIl we summarize the main 
results of the paper. Throughout this work, we take the 
spacetime metric signature to be — \- H — h and we work 
in natural units where c — eq — fiQ ^ I. 



Lorentz violation as dynamical fields inducing sponta- 
neous breaking of Lorentz symmetry, with certain restric- 
tions placed on their dynamics.^ Similar methods can be 
adopted for the matter-gravity couplings as well [20| . The 
linearized equations in this formalism include, as special 
cases, modelsjjf spontaneous Lorentz-symmetry breaking 
with scalar [27J], vector [23], and two-tensor fields |29l.l30l|. 
In linearized gravity the metric is expanded as 



II. FIELD EQUATIONS 



Qfiiy — VfJ^l^ ^~ "'//!>■ 



(3) 



The CPT-even coefficients for Lorentz violation in the 
photon sector of the minimal SME are denoted {kp)^^^'^^ ^ 
which is assumed totally traceless by convention, and 
have all of the tensor symmetries of the Riemann tensor 
and therefore contain 19 independent quantities [9l llllj. 
Following Ref. [13?], it is useful to split these 19 coeffi- 
cients into two independent pieces using the expansion 



{kFY""^^ = C^'"'''^ + -[7^^'''{cfY^ -V^^{cfY^ 



-■q [cfY 



V^'icF) 



(1) 



With this decomposition 9 coefficients are contained in 
the traceless combinations (cf)^" — {^fY'^'^ a ^'^'^ 10 
coefficients are in C^^'^'^^ , which is traceless on any two 
indices. The modified Maxwell equations can then be 
written in the form 






(2) 



where F^^, — d^A 



diyA^ and A^^ is the vector potential. 
This result follows directly from the electromagnetic ac- 
tion of the minimal SME in Minkowski spacetime, when 
the electromagnetic field is coupled in the standard way 
to a conserved four-current j^ = (p, J), and when the co- 
efficients are treated as constants in an observer inertial 
frame. 

In the gravitational sector, the coefficients for Lorentz 
violation are expressed in terms of three independent sets 
of coefficients: t^^'^'^^^ s^'^, u. The t coefficients are taken 
as totally traceless and have the symmetries of the Rie- 
mann curvature tensor, implying 10 independent quan- 
tities. The s coefficients are traceless and contain 9 in- 
dependent quantities. With the scalar u, there are in 
general 20 independent coefficients describing Lorentz vi- 
olation in the gravitational sector. 

Unlike the SME in Minkowski spacetime, it is not 
straightforward to proceed directly from the gravitational 
action to the field equations. This is because introduc- 
ing externally prescribed coefficients for Lorentz viola- 
tion into the action can generally confiict with the funda- 
mental Bianchi identities of pseudo-riemannian geometry 
[lO|. It turns out, however, that spontaneous breaking of 
Lorentz symmetry evades this difficulty [13, [2a| . In Ref. 
[l8| , the linearized gravitational field equations were de- 
rived using a formalism that treats the coefficients for 



Within the minimal SME approach, the linearized field 
equations can be written in terms of the vacuum expec- 
tation values of the coefficients for Lorentz violation, de- 
noted i'^"'" , s'^'', M, which are taken as constants in a 
special observer coordinate system."^ The linearized field 
equations take the form 

-s"^i?^,, + -Sf,^R + Tj^^s^^R^x, (4) 

where Gn is Newton's gravitational constant. In this ex- 
pression Rk^wX is the Riemann curvature tensor, G^^ is 
the Einstein tensor, R^i, is the Ricci tensor, and R is the 
Ricci scalar. All curvature tensors in Q are understood 
as linearized in the fiuctuations hf^^. Since the u coef- 
ficient only scales the left-hand side, it is unobservable 
and is discarded for this work. 

Because of a tensor identity [l^l, the 10 coefficients 

t vanish from the linearized equations, thus leaving 

the 9 coefficients in s'^'^ in this limit. This immediately 
implies that, should an analogy exist between the photon 
and gravity sectors of the SME, it involves a subset of the 
(/cp)^"'^'^ coefficients. This subset is comprised of the 9 
coefficients {cf)^'^- 



III. FIELD MATCH 

A. Conventional GR case 

In GR and Maxwell electrodynamics, the analogy be- 
tween certain components of the metric fluctuations h^^, 
and Afj, reveals itself from the field equations in the har- 
monic gauge: 



d^'hu, = 0. 



(5) 



^ This formalism can be considered a subset of the gravity sector 
of the SME expansion, and was recently dubbed the "Bailey- 
Kostelecky formalism" [2al . 

^ The reader is warned not to confuse the bar notation which in- 
dicates the vacuum expectation values of the tensor fields i^"'"'^ ^ 
s*"^ , and u with the bar notation h^i, used in this paper for the 
trace-reversed metric fluctuations. 



Here h^i, are the usual trace-reversed metric fluctuations 
defined by 



f^^i^ i^^iy nVf^l^'^ c 



(6) 



In the absence of the coefficients for Lorentz violation 
{kp)^^'^^ and s'^'', the Einstein equations in this gauge 
read 



nhu 



-l&TT Gn{Tm)^u, 



(7) 



while the Maxwell equations, in the gauge d^A^ = 0, arc 



aA„, = - 



Jm- 



(8) 



To match the structure of the Maxwell equations one 
typically makes a slow motion assumption for the matter 
source. For example, for perfect fluid matter with ordi- 
nary velocity v^ much less than one, and small pressure. 



(Tm) 00 
(Tm)o] 

{TM)]k 



pv-'v . 



(9) 



Thus, examining Eq. ([7]), it can be seen that the compo- 
nents hjk will be one power of velocity more than hoj, 
and hence negligible. To be more precise, if one adopts 
the standard post-newtonian expansion and counts terms 
in powers of mean velocity v, labeled as 0(1), 0(2), etc., 
one finds from Eq. ^ that 



/iOO '' 


- 0(2), 


hoj - 


- 0(3), 


hjk " 


^ 0(4). 



(10) 

Furthermore, in post-newtonian counting, partial time 
derivatives obey the post-newtonian counting [3l|, 1321 

d V 

dt r ' 

where r is the mean distance. A consistent approxima- 
tion including up to 0(3) terms would take n « V^ and 
Eq. ([7]) would become 

V^/lQ^ = -167rGAr(rA/)o^, 



(11) 



fl2) 



which can be compared with the stationary equations for 
electrostatics and magnetostatics 



V^A, 



-Jm- 



(13) 



From these two expressions it is clear that, given solu- 
tions for Af^ in the stationary limit, the solutions /iq^ 
can be obtained in the manner below. 

1. Replace charge density pq with mass density pm 
and electric current density J^ with mass-current 
density pv^ . 

2. Write down the metric components as /lo/x — 
— IGttGatA^. 

This method agrees with standard results in the litera- 
ture in,!!!. 



B. Lorentz-violating case 

Equations (O and ^ lead to a direct correspondence 
between the solutions for /iqm and A^. In the presence of 
Lorentz violation, this direct analogy involving the trace- 
reversed metric fluctuations disappears because the coef- 
ficients s'^'' in the modified equations ^ generally mix 
the components of ho^ with hjk ■ As a result of this mix- 
ing, hjk contains terms of 0(2) in post-newtonian count- 
ing, in contrast to the GR case ^TU\\ . and so there is no 
particular utility in using the trace-reversed metric fluc- 
tuations hf^i, over the metric fiuctuations /i^,y. 

We focus on the stationary limit, where a match be- 
tween the electromagnetic and gravity sectors can be 
obtained for the metric components /iqo and hoj. This 
gravitoelectromagnetic correspondence is most easily ob- 
tained directly from the stationary solutions to Eqs. ([2]) 
and (jll) for the metric g^i, and the vector potential Ag . 
The gravitational solutions were obtained in Ref. [3 
while the results in electrodynamics were obtained in 
Rcfs. EH, [11]. 

Before displaying the solutions here, it will be conve- 
nient to introduce various potential functions that take 
a similar form for both the electromagnetic and gravi- 
tational sectors. The key source quantities appearing in 
these potentials are the charge (mass) density p and the 
charge (mass) current J-'. The needed potentials are 



U 

lJ3k 
j^jkl 



a 



-4^d^x\ 



p{x'){x-x'y{x~x'Y ,3^, 



?'|3 



a^d^.', 



jHs'){x ~ x')Hx - xy 



In the stationary limit, all partial time derivatives of 
the potentials vanish. The density p is time indepen- 
dent and the current is transverse, djJ^ = 0. This 
implies some simplifications of the identities among the 
potentials listed in Ref. 18], including djV^ = and 

a^-x^'^' = 0. 

The electromagnetic potentials are obtained by inter- 
preting p as charge density, J^ as a steady-state current 
density, and letting the constant a = l/47r. For the grav- 
itational sector, the potentials are obtained by interpret- 
ing p as mass density, J^ — pv^ as mass-current density 
and letting a = Gn- 

The components of the metric fluctuations /iq/^, rele- 
vant for comparison with the electromagnetic sector, can 
be obtained after an appropriate coordinate gauge choice. 
We choose coordinates such that 



djhoj = 0, 



dkhkj — 7:9j{hkk — hoo), 



(15) 



and the metric fluctuations are time independent. To 
post-newtonian 0(3), the metric components Hq^ are 
then given by 



"nn — 



T-Oj 



(2 + 3sO" 



U^gjkjjjk _ 
Ok^jk _ ^n 






(16) 



where a = Gjv is chosen in the expressions (HH). Al- 
though they are not relevant for the match between the 
two sectors, for completeness, the remaining components 
of the metric hj^ are given by 



hjk = [(2-s""')t^ + s""C/""](5- 



0k _^-]l^. 



Ik 



_^^kljjlj 



2sO^>JJ3k 



(17) 



which is valid to post-newtonian 0{2). 

In the electromagnetic sector, we choose the stationary 
limit and adopt the U{1) gauge condition djA^ = 0. The 
modified Maxwell ec^uations have the solutions 



A" = [i+u^FrpE+u^Fyuj,' 



icF^V^ 



A^ 



^C^jok^jk 
'We 

jk-t/k 



f^Ojkl v-ljk 



[CF, 



_/^Okjl Tiki 



E 

f + [i-i{cFr]v^E 

klfvklj 



licFru: 



2^ 
jkh 



/^OjOkjrk 



r^jklm Y^mkl 



(18) 



where the subscript E reminds us to take a = l/An in 
the potentials (HH). 

A glance at Eqs. (fT6|) and ([T8|) reveals that many of 
the same terms occur in both sectors. However, in the 
electromagnetic sector the contributions from the 10 in- 
dependent coefficients C^'^'^^ do not vanish. To match 
the two sectors we must first restrict our attention to the 
special case where 



C' 



ukX 



0. 



(19) 



Next we split the terms appearing in A^^ and /lo/i into 
those involving potentials derived from charge density p 
and those derived from current density J^ . These fields 
are defined as 



ihp)oo - 


= (2 + 3s°")C/ + s^■'=t/^'^ 


{hj)oQ = 




ihj)o3 = 


= -4(1 + is™)V^^' + 2P''V'' 




+2s'^'(X'='^-X^'='), 


(^p)° = 


= [l + UcFr]UE + hi^^Fy'ui,\ 


(Ajf - 


- -{cFrvh 


iA,y = 


- \{cFtWE + \{cFf''V'^^. 


(AjY - 


- [1 - i(c^)™]^i - |(cf)^'=f| 




-\{cEr\xf-x^^\ 



Quantity 


Electromagnetic sector 


Gravitational sector 


Coefficients 


{'^fY" 


gM-- 


Scaling 


\I^-K 


Giv(l + s°°) 


Density p 


charge density 


mass density 


Current J^ 


current density J^ 


mass current J^ = pv^ 


p fields 


(^p)m 


(hp)o^, 


P fields 


kAj)^ 


{h.j)op 



TABLE I: The gravitoelectromagnetic correspondence be- 
tween the electromagnetic and gravitational sectors of the 
minimal SME. 



Note that the split of A^ and h^^ corresponds to splitting 
the terms in the post-newtonian metric into 0{2) and 
0(3) and splitting the terms in the electromagnetic po- 
tentials into "post-coulombian" terms of 0{2) and 0(3) 
[3ll [37| . The correspondence between the two sectors is 
summarized in Table HI 

Given a stationary solution to the modified Maxwell 
equations dH) in the Coulomb gauge {djA' = 0), one can 
obtain the corresponding metric components by using the 
following procedure. 

1. Set Cf"""^ = 0. 

2. Replace {cfY" -^s^''. 

3. Separate A^ into density-sourced and current- 
sourced terms (^p)^ and {Aj)^. 

4. Replace charge density pq with mass density pm 
and electric current density J'-> with mass-current 
density pv'-> . 

5. Write down the metric components 



(Mop = -8^Gjv(l + sO")(^p)^, 
ihj)o^ = -16^G^(l-t-s°")(Aj)^, 

and omit any subleading order terms {0{'s'^)) 



(21) 



(20) 



The close resemblance of the effects of Lorentz viola- 
tion on gravity and electromagnetism is remarkable con- 
sidering the qualitative differences between the theories, 
particularly in the starting lagrangians and field equa- 
tions [i3]- On the other hand, since there is a known 
analogy between A^ and /lop in the conventional case, 
and both sectors are affected by two-tensor coefficients 
for Lorentz violation, one might have expected a close 
correspondence in the appropriate limit. In fact, the map 
constructed above further justifies the construction of the 
post-newtonian metric using the formalism in Ref. jlSf . 
which itself relied on several assumptions concerning the 
dynamics of spontaneous Lorentz-symmetry breaking. 

An interesting feature of the solutions for Lorentz- 
violating electrodynamics is the mixing of electrostatic 
and magnetostatic effects in the stationary limit. As can 
be seen from pp]) , this occurs because a part of the scalar 
potential A^ depends on current density and part of the 



vector potential A depends on charge density, a feature 
absent in the conventional case. This was aptly named 
electromagnetostatics (EMS) in Ref. [33] ■ For Lorentz- 
violating gravity, a similar mixing occurs and /iqo de- 
pends partly on mass current while Hqj depends partly 
on mass density, resulting in what can be called gravito- 
electromagnetostatics (GEMS). These features are illus- 
trated with specific examples in Sec. [V] 

Note that other possibilities are open for exploration 
concerning the match between the two sectors of the 
SME. For example, we do not treat here the interesting 
possibility of whether an analogy persists using gravita- 
tional and electromagnetic tidal tensors, as occurs in the 
Einstein and Maxwell theories 13811. 



IV. TEST-BODY MOTION 

In this section we study another aspect of gravitoelec- 
tromagnetism. This concerns the behavior of matter in 
the presence of the stationary gravitational or electric 
and magnetic fields. As we show below, if one adopts the 
appropriate limit, the behavior of test masses in gravi- 
tational fields and test charges in electric and magnetic 
fields is analogous, despite the presence of Lorentz vio- 
lation. However, differences do arise in the presence of 
Lorentz violation when comparing the gravitational spin 
precession to the classical spin precession of a magnetic 
moment in the presence of electromagnetic fields. 



A. Geodesic motion 

When Lorentz violation is present in the electromag- 
netic sector only, test charges e obey 






m 



(22) 



where u^^ is the four- velocity. With the usual identifica- 
tion of the electric and magnetic fields, Ej — Fjq and 
Bj = {l/2)ejkiFki, we can write the spatial components 
of ([2^ as the familiar Lorentz-force law: 



^ ^ ^[E^ + fiJ y, BV] 
at m 



(23) 



For small velocities, u^ Ki v^ — dx^ /dt. Thus, with kp 
affecting only the electromagnetic sector, the force law 
for charges is conventional jll| . 

In the SME, restricted to only the s^'' coefficients, 
freely falling test bodies satisfy the usual geodesic equa- 
tion 



dui^ 
~d^ 



'^\p^"^' 



(24) 



In its full generality, the structure of (|24p is quite different 
from Eq. (1^^ for charges. Nonetheless, in the weak-field 
slow motion limit of gravity, there is a correspondence. 



Changing variables in ((24)) to coordinate time, one can 
solve for the coordinate acceleration a^ = dv^ /dt in terms 
of the connection coefficients projected into space and 
time components using standard methods. One obtains 
the well-known expression [32l |. 



-(r%o + 2rV' + r°,,«S>^-. (25) 



So far, Eq. (P5|) is an exact result, and bears little 
resemblance to Eq. (l23l) . If one then assumes that the 
test particle velocity is small and keeps only terms linear 
in the test particle velocity v^ , the acceleration becomes 



^ 00 ^^ Ok" 



J- 00^ ■ 



(26) 



To get a match with equation (|23p additional assump- 
tions are needed. For example, in the post-newtonian 
approximation, the dominant contributions to the con- 
nection coefficients are given by the formulas 



-L n 



dogok 
1 



1 



g^'^dkOoo, 
1 



Ok 



-L n 



9o9jk + -^{dkgoj - djgok), 



1 



2^0500, 



(27) 



which is valid to post-newtonian 0(4). If the metric is 
stationary in the chosen coordinate system, (i9o5/ii/ = 0), 
the acceleration, in terms of the metric fluctuations 
hf,i, = 9^,1, - rjf^^, is given by 



-djhoo + v^(djhok - dkh^j) 



1. 



-hjkdkhoo, 



(28) 



which neglects terms proportional to the test-mass ve- 
locity squared but otherwise is valid to post-newtonian 
0(4). This expression now resembles the Lorentz-force 
law ([23l) except for the last nonlinear term. 

To be consistent with the post-newtonian approxima- 
tion to 0(4), the last term must be included, as well 
as nonlinear contributions to /iqo at 0(4). This is be- 
cause the second term in Eq. ([^5]). the so-called grav- 
itomagnetic acceleration term, is an 0(4) term in the 
post-newtonian expansion. 



1. GR case 

Results from GR are contained in (PSI) and ((231) in the 
limit of vanishing coefficients for Lorentz violation. In 
the stationary limit of GR, and in the coordinate gauge 
(flSl). the acceleration (l28l) can be written as 



dj4>~Av''{djV'' -dkV^). 



(29) 



Here is a post-newtonian potential that includes 0(4) 
terms in GR IH: 



{p + pU + 3p-2pU)^,^, „ 2^.^ 



(30) 



where p is the perfect fluid pressure and 11 is the internal 
energy per unit mass. Note that (p does not satisfy the 
field equation (fT2|) . 



vhoo = —IGttGnP- 



(31) 



Instead it satisfies 



V^^ = -47rGjv ip + pn + 3p- 2pU) - 4( VC/)^ (32) 

Therefore (f> 7^ ^ooi and it cannot be obtained directly 
from the solutions to Aq in Eq. ([T3|) using the standard 
match. 

Generally, care is required in discarding the nonlinear 
terms in (/), while keeping the second, gravitomagnetic 
terms in Eq. (l29l) . A simple estimate for a realistic sce- 
nario can establish this. For a rotating spherical body, 
the solution for V^ is of order GnI^/t'^ ^ GNMR^uj/r"^, 
where / is the inertia of the body, R its radius, uj its an- 
gular velocity, and r is the coordinate distance from the 
origin to the location of the test body. The typ ical test 
particle velocity v^ is of order v ^ ^Gj^Mjr or less, 
where M is the mass of the source body. Thus, the con- 
tribution to (j29p from the gravitomagnetic force term on 
a test particle outside of the source body, has an approx- 
imate size 



{GnM)R'^ujv 



^gm I 



(33) 



The contribution from the nonlinear terms in (j) to the 
test particle acceleration have an approximate size US/U 



GmMGmM 



\ani\ 



(34) 



Assuming that the nonlinear contributions are much 
smaller than the gravitomagnetic contributions, |a„i| << 
|agm|, amounts to assuming 



Rujv » 



GnM 
R 



(35) 



For example, consider a test-body near the Earth's sur- 
face. For this case one finds that condition (|55|) implies 
the unrealistic condition that the test particle velocity 
must be greater than 1/2000 of the speed of light. 

In addition to the above argument, it is important to 
recall that terms of second and higher order in the test 
body velocity v^ were discarded in (|26p . In terms of post- 
newtonian counting, these terms make contributions to 
the acceleration a^ at the same order (0(4)) as the non- 
linear terms. One example is the term F-'^.^ii'^f', which 
can be shown to have an approximate size similar to p4p 



in the typical post-newtonian scenario.^ Furthermore, it 
is interesting to note that an argument along the lines of 
the one presented here appeared in the original paper by 
Lense and Thirring in 1918 2]. There it was emphasized 
that nonlinear terms must be included in the equations of 
motion, in addition to the gravitomagnetic force terms, 
to properly account, for example, for the precession of 
the orbital elements of the planets. As an alternative to 
this reasoning, one can incorporate the nonlinear terms, 
such as those occuring in Eq. (P^ , to form "Maxwell- 
like" equations, as pursued in Ref. (39| . 

For simplicity here we separate out the gravitomag- 
netic and gravitoelectric acceleration terms from the non- 
linear terms. Thus we write 



a — acEM + a-NL, 
where the separate terms are given by 

Sgem = Eg + v'xBg, 



(36) 



(37) 



Here we have identified the gravitoelectric and gravito- 
magnetic fields for GR: 



Eg - Vt/, 

Bg ^ -4V X V. 

2. Lorentz-violating case 



(38) 



To see if there is any resemblance for the Lorentz- 
violating case between the gravitational force law and the 
electromagnetic force law, we can proceed from Eq. (j26p . 
Adopting the stationary limit (1^51) . we restrict attention 
to the gravitoelectromagnetic portion of the acceleration 
which we denote (oOgba/- This acceleration is given by 



{a'y 



GEM 



-djhoo + v''{djhok - dkhoj), (39) 



For a consistent expansion to first order in the coefficients 
s^" , we take hoo to 0(3) and Hqj to 0(2). This produces 
an acceleration to first order in the coefficients s'^^ that 
is at most 0(3). 

In the presence of the coefficients for Lorentz violation 
s^"^, the components of the metric from Sec. IIII Bl are 
needed to this order: 



'loo — {hp)ao + (^.7)00: 
ho J = ihp)oj- 



(40) 



Note that the expansion of hoj is truncated at 0(2) since 
this term is multiplied by a velocity (0(1)) and therefore 
produces an 0(3) term in the acceleration. 



^ In fact, in the case of laboratory gravitational sources and test 
bodies, these velocity-squared terms may be substantially larger 
than the nonlinear terms in Eq. I I29II . as discussed in Ref. [4C|. 



With the considerations above, the gravitoelectromag- 
netic acceleration can be written to 0(3) as 



{a')GEM = Eg + v X Bg, 



(41) 



which now resembles the result in Eq. (P5|) . The effective 
electric and magnetic fields are given by 



1 






(42) 



This result demonstrates that in the limit that the gravi- 
toelectromagnetic acceleration terms are considered, the 
force on a test body takes the same form in the electro- 
magnetic and gravitational sectors of the SME. 

To use this result in a manner consistent with the post- 
newtonian expansion, additional terms at 0(4) but at 
zeroth order in the coefhcients s'^'^ need to be included 
in the acceleration. Specifically, the total acceleration at 
0(4) takes the form 

a^ = {a'y + a^NL + «'[5,(/ij)ofc - dk{h.j)aj], (43) 

where onl is given by Eq. (1371) and the components 
{hj)oj are taken to zeroth order in the coefficients s^"^. In 
the limit s'^'^ = 0, this expression reduces to the standard 
GRresuh in (l36l). 



B. Spin precession 

The classical relativistic behavior of a particle with 
a magnetic moment jl under the influence of external 
electric and magnetic fields is well known. Consider a 
particle, such as an electron, with spin s defined by 



Ai 



2m 



(44) 



Here, e is the charge of the particle, m is the mass, and 
g is the gyromagnetic ratio for the particle. We can 
describe the behavior of the spin relativistically using 
the spin (spacelike) four- vector S^ which, in an instanta- 
neous comoving rest frame, takes the form (S"" = 0, S^ = 
s^ ) . The motion of the particle is described with the four 
velocity u^, which satisfies Eq. (1^^ . In addition, we have 
the identity S^^Ufj, = 0. 

If we ignore field gradient forces and nonelectromag- 
netic forces, the behavior of the classical spin four-vector 
S^ is determined by the dynamical equations [J, |4l| 






e 
m 



^F'^^'S^ 



^-l)u^iS,F''^ux) 



(45) 



A formula for the precession of the spin as measured in 
a locally comoving reference frame can be obtained by 
projecting S'^ along comoving spatial basis vectors e^-, 

and making use of Eqs. (H5|) and (1^^ . With the choice of 



(7 S3 2, the lowest order contributions to this precession 
can be written 



dr 



1 



S X B - -S X {v X E) 



k 



kj- 



(46) 



This result holds up to order u^ in the particle's ordinary 
velocity. Furthermore, Eqs. (|^5|) and (|46p will still hold 
in the presence of Lorentz violation in the photon sector 
since the force law takes the conventional form ([23]) . 

The behavior of the classical spin four-vector in the 
presence of gravitational fields is given by the Fermi- 
Walker transport equation [42] 






= -r'^,X^ +"''K^-)' 



(47) 



where a'^ is the acceleration of the spinning body. For 
comparison with the electromagnetic case, we assume 
that the spin is in free fall {a'^ = 0), and again find 
the spin precession along the comoving spatial basis e^-, 

a standard technique [3l|, \^ . The resulting precession 
was obtained in the post-newtonian limit for an arbitrary 
metric in Ref. [18] and is given by 

-^ = ^kjS" [kiv'^djhoo - v^dkhoo) + hid.hok - 9,.^) 

+ ^v'{djhki-dkhji)]. (48) 

which is valid to post-newtonian order 0(3). Since this 
result was derived for an arbitrary post-newtonian met- 
ric, it holds for the metric in Eqs. (fT6| and (fT7|) as 
well. Note that the expression (j48|) does not immedi- 
ately match (|46p due to the last terms in (P5| dependent 
on hki at 0(2). However, a judicious choice of coordinate 
gauge may alleviate the problem, as we show below. 

In GR, we can make use of the results of Sec. IIII Al in 
the harmonic gauge. When expressed in terms of h^i, the 
GR spin precession to 0(3) is 



d_S^ 
dr 



S- X (V X .g)*^ - Sx{vx V/iQo)* 



^kj, 



(49) 



where g^ = hoj and we have omitted contributions from 
hjk ^ 0(4). The expression (|49|) now resembles the elec- 
tromagnetic counterpart, at least up to numerical factors. 
In fact, one can again define effective electric and mag- 
netic fields for gravity: Eg = (l/4)V/ioo, Bg = V x g. 

We next introduce Lorentz violation in the gravita- 
tional sector in the form of the post-newtonian metric 
(fTHl) and ([TT]) . Unlike in GR there are off-diagonal terms 
in hjk that cannot be eliminated by a choice of coordi- 
nate gauge. As a result, we find that the third term in 
(|48l) cannot be reduced to a term of the form S x V$, 
where $ is a scalar. Therefore it is not possible to match 
the form of the spin precession in the gravitational sec- 
tor to the electromagnetic sector of the SME, the latter 
of which takes the form (j46|) . Evidently, this is due to 
the important role of the metric components hjk in the 
general spin precession expression (|48p . 
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V. EXAMPLES AND APPLICATIONS 

In this section we illustrate the methods of matching 
electromagnetic solutions for the fields to gravitational 
solutions for the metric components. We also demon- 
strate the match between the two sectors for test-body 
motion. In our examples we study both a static pointlike 
source and a rotating sphere. Finally, we comment on 
the observability of the GEMS mixing effects in specific 
gravitational tests. 



A. Static point source 

We consider first a point charge q at rest at the origin 
in the chosen coordinate system. The potentials in the 
Coulomb gauge were obtained in Ref. [35[ and are given 
by 



{kF)°^°^ - {kpf^^'^x^x''] 



A° = ^- [1 

Anr '- 



(50) 



where x = x/r and r = \x\. 

Using the method outlined in Sec. IIIIB| we can obtain 
the corresponding metric components ho^ in the fixed 
coordinate gauge (fT5|). First we expand the coefficients 
{kpY^^^^ into C and cp terms using ([1]). Next, we set 
all of the coefficients C = 0, according to step I. Then 
we make the replacement in the remaining coefficients 
Cp — > s. At this intermediate stage the potentials are 
given by 



A" = -^ [1 
47rr '- 

q_ 

8T:r 



I5OO + ^s^'^x^x''] 



Ai = -^[^i+^^x'^x'] 



An interesting feature arises from this simple solu- 
tion. In Lorentz-violating electromagnetism, even a 
static source will generate a magnetic field. For gravity, 
the analog of this effect occurs. For example, consider 
the scenario in which the coefficients s^*^ = 0. Apart 
from a scaling, the gravitoelectric force appears conven- 
tional. However, even when the source body is static, a 
test body with some initial velocity wq will experience a 
gravitomagnetic force. 

The nature of this gravitomagnetic force is illustrated 
in Fig. [T] The gravitomagnetic field itself falls off as 
the inverse square of the distance from the point mass, 
and curls around the direction of the vector denoted s, 
where s^ = — s"-'. A test mass approaching the pointlike 
source will be deflected in the opposite direction of s, as 
illustrated in the figure. 



i 



4 ^ 



o 

t 






FIG. 1: The gravitomagnetic field Ba (blue arrows) from 
a static point mass m (center). The field curls around the 
direction of s (red arrows) and falls off as the inverse square of 
the distance. An approaching test body is deflected opposite 



(51) 



Since there is no dependence of the potentials on any 
current density, for step 3 we simply note that in Eq. 
^ A^ = [Apf and A^ = [ApY . We make the re- 
placement q ^ m and multiply the potentials by a factor 
of — 87rGjv(l + s"") and cancel subleading order terms 
{0{s'^)). This yields 



2GNm r 

hoo = 1 

r '■ 

Gntti r_^ 



1^00 + Is^'^x^x''] 



T-Oj 



:Oj+sOfe^fc^j 



(52) 



In a similar manner, we can also obtain effective grav- 
itoelectric and gravitomagnetic fields using (|42l) : 



^_^ ^ _^[^.(1 






r 
26* Arm 



3-00 
2* 



'{E^^'x^x') 



e^'^'s°'=x'. 



(53) 



Using these expressions the acceleration of a test mass 
can be written in the Lorentz-force law form (|4ip . 



B. Rotating sphere 

We next turn our attention to a more involved exam- 
ple, a spherical distribution of charge or mass that is 
rotating. In Ref. [35|, a scenario was considered that in- 
volved a magnetized sphere with radius a and uniform 
magnetization M. In conventional magnetostatics, an 
idealized scenario would allow for the sphere to have zero 
charge density and no electrostatic field surrounding it, 
thus it would only produce a dipole magnetic field. In 
the presence of Lorentz violation, however, a dipole elec- 
tric field persists, with an effective dipole moment con- 
trolled by the parity-odd coefficients for Lorentz violation 

{kpf^'^K 

Since we aim to find the gravitational analog of this 
solution, we cannot consider an object with zero charge 
density. Instead we study a closely related example: 
a charged rotating sphere, which produces an effective 
magnetic dipole moment rfi in the conventional case. For 
this example, the current-induced portion of the electric 



scalar potential, {Aj)°, can be obtained directly from Eq. 
(31) in Ref. [H: 



^ £j'='(cF)°^xfem' 

47^7-2 



(54) 



which holds for the region outside the sphere. For a ro- 
tating charged sphere 



m^ = i/fiW-' . 



(55) 



where tS is the angular velocity of the sphere. The quan- 
tity Ie is the charge analog of the spherical moment of 
inertia for massive body, 



Ie^ d^xp\x\ 



(56) 



Comparing (|54p with the standard dipole potential, the 
effective dipole moment is 



p^ = f^^^m^icF) 



01 



(57) 



The effective electric field therefore takes the standard 
form 



E = 



3p ■ XX — p 
47rr^ 



(58) 



The gravitational analog for the solutions (IMl) and ([55)1 
can be obtained using the methods in Sec. IIIIBl Since 
the C'^^^'^ coefhcients do not appear, step 1 is redun- 
dant. We next make the replacement [cpY'-' — > '^^ ■ All 
that remains is to change Pq -^ p„i and multiply (1541) by 
IQttGn which yields 



[h 



./J 00 



AGNle^'''x^uj''s°^ 
3^ ' 



(59) 



where now I is the spherical moment of inertia of the 
massive body, given by Eq. (j56p using mass density. Note 
that this produces an extra component of the gravitoelec- 
trie field Eg = {l/2)V{hj)oo- 

In the electromagnetic case, part of the electrostatic 
field arises from the effective current of the rotating 
charged sphere, a feature absent in the standard Maxwell 
theory. This unconventional mixing of electrostatics and 
magnetostatics has an analogy for stationary gravita- 
tional fields produced by a rotating mass, in the presence 
of Lorentz violation. Thus a uniformly rotating sphere of 
mass produces a gravitoelectric field whose strength de- 
pends on the rotation rate, a feature absent in standard 
GR. 

As in the point- mass example, the vector s is responsi- 
ble for the effect. In Fig. [51 the effective dipole moment of 
a rotating spherical mass is depicted. The dipole moment 
is obtained from the cross product of s with 4/^/3. 

The full solution for the case of a rotating massive or 
charged sphere can be constructed using the potentials U, 
U^'^, V^ , and X^^'- in Eqs. (ITl)) . For the electromagnetic 



i 

t 



FIG. 2: A depiction of the effective dipole moment that de- 
velops for a rotating sphere in the presence of the coefficients 
for Lorentz violation s. The dipole moment (green arrow) is 
proportional to the cross product of s (red arrows) with the 
angular momentum of the sphere (blue arrow). 



case (a = l/47r) we obtain, for the region outside the 
sphere r > R, 



Ue 



u- 



jk _ 



V^ 

^E 



Q 

47^7-2 ' 

Qx^x'^ 



h 



{5^^ - ix^x^) 



127rr2 



X 



jki 



W^ 



-k -J 



I'e 



/^J*^' 



/£;(eJ'="x'w" + eJ'™x'=w™) 
127rr2 



31' 



Met' 



(60) 



where I'^ is a spherical moment given by the integral 



in Eq. ([55]) with \x\ instead of 



Using these ex- 



pressions it is straightforward to calculate the associated 
electric and magnetic fields as well as the gravitoelectric 
and magnetic fields. The expressions are lengthy and 
omitted here. 



C. Applications 

A full analysis of the dominant observable effects in 
gravitational experiments and observations has been per- 
formed in Ref. 118!]. However, the coefficients were ana- 
lyzed collectively and the separation of various distinct 
Lorentz-violating effects was not fully studied. Here we 
focus specifically on the observability of the novel gravit- 
omagnetic force shown to arise in the point-mass example 
in Sec. IV A] and illuminate its role in a key test. 

Lunar laser ranging and atom interferometry have 
measured 8 of the 9 coefficients in s'^'^ and the combined 
results are tabulated in Ref. [21|. These results are re- 
ported in the standard Sun-centered celestial-equatorial 
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frame (SCF), where coordinates are denoted with capital 
letters for clarity. In this frame, the current constraints 
on 's'^^ are at the 10~^ level. For 'S^^ ^ the constraints 
are at the weaker level of 10^^-10^^. The gravitomag- 
netic force due to the effective gravitomagnetic field in 
the second of Eqs. (|53l) is controlled by the 's'^'' coeffi- 
cients. This force has been measured by both lunar laser 
ranging and, effectively, atom interferometry. However, 
its specific effects are most easily discernable in orbital 
tests such as the lunar laser ranging scenario, so we focus 
on this case. 

The principle effects from the s coefficients for lunar 
laser ranging are modifications to the relative accelera- 
tion of the Earth and Moon. This acceleration includes 
such terms as the gravitomagnetic terms considered in 
Eqs. ((55)) . In fact, from the results in Ref. [ISj], one can 
read off the portion of the Earth-Moon acceleration 5a^ 
responsible for the effective force that is described in Fig. 
[TJ In the SCF coordinates, it reads 



J 2GAr(5m K/-TK J 

da = ^ V [s r 



5^'V^) 



(61) 



where 6m is the mass difference between the Earth and 
Moon, r"^ is the coordinate difference between the Earth 
and Moon center of mass positions and v' is their relative 
coordinate velocity. 

The dominant observable effects from Eq. ((6T|) are os- 
cillations in the lunar range at the mean lunar orbital 
freqeuncy lo. In the lunar laser ranging scenario, these os- 
cillations are controlled by two linear combinations of the 
s^"^ coefficients called s°^ and s"^, which are expressed 
in the mean orbital plane of the lunar orbit. These two 
quantities control the size of the Lorentz- violating grav- 
itomagnetic force for this case. Using over three decades 
of lunar laser ranging data, analysis reveals that s"^ = 
(-0.8 ± 1.1) X 10-6 and s°^ = (-5.2 ± 4.8) x 10"^ ^. 
Therefore there is no compelling evidence for the gravit- 
omagnetic force controlled by s'^"' coefficients. However, 
ongoing tests such as the Apache Point Observatory Lu- 
nar Laser-Ranging Operation have already improved on 
lunar ranging capability and could significantly improve 
sensitivity to this effect (43,]. 



VI. SUMMARY 



In this work we have shown that an analogy exists be- 
tween the gravitational sector and the electromagnetic 
sector of the SME at two levels. First we showed that 
in the stationary limit and for a particular coordinate 
choice, part of the post-newtonian metric /iq^i in the grav- 
ity sector can be obtained from the vector potential A^ 
in the electromagnetic sector by essentially making a se- 
ries of substitutions, most notably the exchange of the 
coefficients C^"" -^ s^"" , as outlined in Sec. lHIBl For the 
equations of motion of a test body, the gravitational case 
was shown to resemble the electromagnetic Lorentz-force 
law, so long as nonlinear terms in the geodesic equation 
are disregarded. 

In Sec. |Vl we provided two examples of how the mixing 
of electrostatics and magnetostatics in Lorentz- violating 
electrodynamics has an analog in the gravitational case. 
In the same manner as a point charge produces a mag- 
netic field in the presence of the electromagnetic coeffi- 
cients C°^ , we showed that a point mass will produce a 
gravitomagnetic field controlled by the coefficients s''-'. 
Similarly, we also explored the converse of this exam- 
ple, demonstrating that a moving mass produces an ad- 
ditional gravitoelectric field. We also discussed the ob- 
servability of the gravitomagnetic force controlled by the 
s'^-' coefficients in lunar laser ranging tests. 

Several areas are open for future investigation. One 
possibility is to systematically isolate the GEMS mixing 
effects from others in the various predicted signals for 
Lorentz violation in gravitational experiments [18| , along 
the lines of the discussion in Sec. IV CI It also would be 
interesting to investigate whether any analogy is possi- 
ble in the presence of the matter sector coefficients that 
play a role in gravitational experiments [2a, |4J]. Fur- 
thermore, using a method similar to the one developed 
in this paper, it may be possible to extend the class of 
signals for Lorentz violation by looking for gravitational 
analogs of the nonminimal electromagnetic sector of the 
SME, which goes beyond the minimal {kp)'^^^^ coeffi- 
cients [isli. 
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